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Abstract. We discuss the dynamic behavior of a tagged particle close to a 
classical localization transition in the framework of the mode-coupling theory of 
the glass transition. Asymptotic results are derived for the order parameter as well 
as the dynamic correlation functions and the mean-squared displacement close to 
the transition. The influence of an infrared cutoff is discussed. 
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1. Introduction 

Recently, there has been renewed interest in the slow transport of particles in 
heterogeneous disordered environments, for example, the diffusion in porous media 
[IHS], or anomalous-diffusion phenomena in binary glass forming mixtures [SHS]. 
Applications arise in many diverse subjects of science and engineering, notably the 
physics of transport in crowded biological systems |9, , ion transport in glass formers 
PUl [TT] , or the geophysics of volcano eruptions [H] . 

One way of modeling such systems is by binary systems composed of a slowly 
relaxing "matrix" component and a "fast" species responsible for the transport. This 
naturally leads to binary mixtures of particles with disparate sizes, be they soft [Sj 
or hard spheres |13l 114] . If the slow matrix is driven to glassy arrest, it forms a 
random heterogeneous background that can be considered as frozen on the time scale 
of fast-particle transport. One can also freeze the matrix from the outset, and is led to 
quenched-annealed mixtures [T5l - flT] or, abstracting even more, the Lorentz gas (LG) 
model [3H5] , a single mobile particle in a random disordered matrix. 

All these model systems show intriguing transport phenomena when the matrix 
is dense enough, although the precise connection for example between the LG model 
and glass-forming binary mixtures remains to be established [18]. Most notably, the 
mobile species can undergo a localization transition, where its long-range transport 
ceases because the embedding matrix is too dense. For the LG, it is understood that 
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this transition is a dynamic critical phenomenon il9) associated with the percolation 
of void space between the frozen obstacles. This leads to distinct scaling predictions 
close to the transition, some of which have only been worked out recently [20] ■ 

One of these predictions is anomalous power-law diffusion: the mean-squared 
displacement (MSD) of the mobile particle, (5r^(i), for t -> oo does not grow linearly 
in time as expected for ordinary diffusion. Instead, 5r^{t — > oo) ~ i^/^ with a dynamic 
critical exponent 2; > 2 is predicted and verified in extensive computer simulations for 
the LG [3]. In three dimensions, z = 6.25 is established, so that 6r^{t — > 00) ^ t^-^^ . 
It is not yet clear, whether this asymptotic prediction also holds for a binary mixture, 
where correlations and residual thermal motion of the background matrix, and the 
excluded- volume interaction of many tracers come into play [6]. 

For dense, glass-forming binary mixtures on the other hand, the mode-coupling 
theory of the glass transition (MCT) provides an excellent theoretical framework 
that allows many detailed predictions [T51 [HI The glass transition in this 

framework comes about because of collective local nearest-neighbor caging; quite 
unlike the single-particle dynamic critical phenomenon that drives the localization 
transition in the LG model involving a divergent length scale. Gonnected to this, the 
glass transition is a discontinuous dynamic transition - the long-time limit (called 
nonergodicity factor) of the collective dynamic density correlation functions jumps 
from zero to a nonzero value at the glass transition. In contrast, the localization 
transition is continuous in this sense and it belongs to a universality class of its own 
[23] . Theories along the lines of the current MGT have been applied to the LG 
previously [24H26] . and the appearance of continuous transitions in MGT applied to 
both quenched- annealed [T5HT7] and size-disparate mixtures [131 [H] cast new attention 
to the description of such divergent-length localization transitions within that theory. 

Most previous studies of localization transitions within MGT have put the 
influence of a divergent length scale on the theory's vertex out of focus, explicitly 
dropping dependencies on wave number gr - in the memory kernel only [24j , or fully in 
the form of a schematic model - or by silently introducing small-g cutoffs through 
numerical discretization. In fact, such a cutoff can be argued for since the small-g 
structure of MGT's vertex is incorrect at least in the LG [55]. However, it seems to be 
largely undocumented so far how sensitive the expected results are to such a cutoff. 

Here we document the dependence of asymptotic results of MGT close to 
localization transitions on the small-wave-number regime. The analytic results are 
corroborated by careful numerical evaluation of the MCT integrals for the tagged- 
particle correlators. We identify the regime where the solutions are dominated by the 
small-g features by studying the dependence of the solutions on an ad-hoc infrared 
cutoff. For simplicity, we obtain numerical solutions for the LG model, but the generic 
picture is the same in the MGT for mixtures or quenched-annealed systems. 

2. Methods 

Consider a mixture of classical hard-sphere particles in d dimensions, with diameters 
da, where Greek indices label the different species. N — is the total number 

of particles in a box of volume V , by Xa = N^/V we denote the respective number 
concentrations. Since we are interested in single-particle transport quantities, we 
consider the dynamic tagged-particle density correlation function, 

Mq,t)^{Qi*{q,t)gi{q)) , (1) 
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where t) — e-xp[iq-fa,s{t)] is the corresponding tagged-particle density fluctuation 
to wave vector q, and fa.s{t) the trajectory of the singled-out tracer particle (taken 
to be of species a). 4>a{q,t) is real- valued and does not depend on the direction 
of the wave vector q for a translational-invariant and isotropic system; we also 
omit the t = Q argument since we assume all dynamics to be time-translational 
invariant. For convenience, it serves to introduce the Laplace transformed 4>a{q, z) = 
i e'^^*'(j)a[q,t) dt for complex frequency z such that the integral converges. 

Starting from Newton's equations of motion, the time evolution of the tagged- 
particle correlation function can be represented by a Mori-Zwanzig equation |21j . 

-1 



ba{q,z) = - 



^hth.a + + rriaiq, z) 



(2) 



where i'th,Q = \/kBT /rria is the thermal velocity of the particles with 
mass mjj, characterizing the ballistic short-time motion. The term Va — 
p X/^gaisRafs \/'^fJ'ai3 ariscs because for hard spheres, the dynamics is 
described by a non-self-adjoint pseudo-Liouville operator , and is a central quantity 
calculated in Enskog theory [30, 31 . For brevity, we have set Rap = da + dp and 
liap — ITT- p/i'm-a +1^1/3), and denoted by gap the value of the radial distribution function 
for two particles of species a and /? at contact. In the time domain, Eq. ^ reads 

^ d^(t)a{q,t)+q^I^adt(l)a{q,t)+q^<l)a{q,t) (3) 



2 '^t 
"th,a 



niaiq, t — t')dt'4'a{q, t') dt' = , 

with initial conditions (j)a{q,0) = 1 and <j)a{q,0) — 0. We are not interested in 
the short-time dynamics. Therefore, numerical solutions are later shown for the 
overdamped case, where the second derivative in Eq. ^ and thus inertia terms are 
dropped. We can then set z/q = 1 to define the unit of time. For the resulting Brownian 
dynamics, the low-density asymptote of the LG has been worked out exactly [32]. 

The memory kernel ma{q,t) is formally given as a correlation function of 
generalized fluctuating forces. MCT assumes the dominant contribution in the 
relaxation of these forces to be given by their overlap with density-pair fluctuations 
- for the tagged-particle dynamics, with pair modes of the form g''^{k,t)g~f{p,t), 

formed with the collective number-density fluctuations ga{q, t) = X^fc^i ^^P['^q''^a,k{t)]^ 
where fa^k{t) is the trajectory of the k-th particle of species a. The corresponding 
collective partial number density correlation functions are denoted by ^apiq,t) = 
{Q*a[q^'t)Qp{q))^ and their initial value is the static structure factor, S{q) = $(q, 0). 
The MCT expression for the tagged-particle memory kernel then reads 

/d'^k 
j—j^ {eqkfcaa' {k)Caa" {k)^ a' a" (k, t)(j)a (p, t) . (4) 

where g is the total number density of the d-dimensional system. It is hence expressed 
as a bilinear functional of the tagged-particle density correlation function 0^ (p, t) , 
and of the collective density correlators ^apik,t). Momentum conservation implies 
q = k+p, and we denote by e^ the unit vector in g-direction. The coupling coefficients 
are given in terms of the Ornstein-Zernike direct correlation functions Cap (q) ■ 

To determine the full dynamics of the system, one would first need to calculate 
the collective dynamics, t), for which analogous MCT expressions can be derived. 
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In this paper, we are only concerned with the tagged-particle dynamics, in the case 
where the slow relaxation time of the heterogeneous matrix, r^^, is much larger than 
the time scales characterizing the nontrivial dynamics of the tracer. We are also 
interested in slow dynamics beyond the initial short-time transient, i.e., in times 
l/i'a ^ t ^ Tap- In this case, ^{q,t) F{q) >~ 0, a positive-definite non- vanishing 
matrix of nonergodicity factors [33]. It is standard MCT procedure to calculate F{q) 
from the given static structure factor of the system. Generically, F{q) decays at large 
q to provide an implicit large-g cutoff of the integrals; it is also non-singular at all q. 



^(^) = Cat3{k)ca^{k)Fp^{k) > (5) 
a/37 

is the effective coupling coefficient that is positive and non-singular. The expression 
of the MCT memory kernel then reads 

/d'^k 
j^^{e^kfv{k)<j)aip,t). (6) 

Rewriting the integral in terms of bispherical coordinates, i.e., in terms of moduli k 
and p, one realizes that the /c-integral can be carried out immediately, and one is left 
with a linear functional 



ma{q,t) = m[(l)]{q,t) = g J dpVqp4>a{p,t) , (7) 

where we omit the subscript labeling the species for brevity. 

The theory given by Eqs. ([3| and (|6| is a special case of MCT for the dynamics of 
a tagged particle in a mixture. In the presence of quenched rather than self-generated 
disorder, the derivation of such a theory becomes more subtle [Mj. Nevertheless, the 
resulting tagged-particle equations are structurally equivalent 17 , as is the MCT for 
the LG, and also the limit of MCT considering the remaining dynamics of particles 
deep inside the glass, the strong-coupling limit [53]. To simplify the further numerical 
procedure, we consider the LG model from now on, i.e., a point tracer among randomly 
placed (overlapping) hard-sphere obstacles with radius R. The vertex v(k) can then 
be evaluated explicitly; in d = 3 one obtains 

vi^oik) = ^{2^fR\3i{kR)f > 0, (8) 

with j'l the first modified spherical Bessel function. Also, v = y/ngR^vth in this case. 

The nonergodicity parameter f{q) = limt_j.oo 0(9,0 can be determined |33ll35j as 
the largest positive solution of 

^^^=m[/](<z)/<z^ (9) 

It is well known that bifurcations exist in these equations: smooth changes of 
the coupling coefficients (such as induced by a smooth increase in density without 
structural changes) can lead to singular changes in f{q). Standard MCT considers 
m to be a nonlinear functional of the /, giving rise to so-called Ai bifurcations (as 
classified by Arnol'd) that one identifies as ideal glass transitions where m{q) jumps to 
a finite value m'^{q) > [ST]. In the present context, we are concerned with continuous 
transitions, viz. the critical points in parameter space, where the trivial solution, 
f{q) — for all finite q, changes continuously to a non-trivial solutions f{q) > 0. 
At the critical point, /^(g) = for all finite q, and consequently m'^{q) = there. 
Note that for q = 0, /(O) = 1 holds always. These continuous transitions arise in 
the LG upon changing the scatterer density p, and in binary-mixture systems through 
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tuning the parameters such that the combination of direct correlation functions and 
nonergodicity parameters entering v{k) leads to sufficiently weak coupling. 

Not all continuous MCT transitions are of the localization-transition type 
discussed here. For example, in dumbbell-shaped molecules with top-down symmetry, 
a glass may form where the rotational degrees of freedom associated with this 
symmetry remain ergodic, and eventually freeze by further compression f36H41| . This 
MCT transition is continuous but does not involve a divergent static length scale. 

If m[f] is a linear functional as in Eq. 0, 

we can express / via Eq. ([9| to get 

m{p) 

+ m{p) 

a closed nonlinear integral equation for m{q), where we have introduced Wqp = Vqp/p^. 
Expressing m in this form, one recognizes immediately that, whenever m'^ is zero 
for finite fc, the peculiarities of the integral at small k become important since 
there is a dangerous denominator. In the following, we will present both analytical 
considerations valid asymptotically close to the transition, and numerical results. The 
latter have been obtained by a specially devised numerical grid for the discretization 



m{q) ^ g I dpVqp ^ ^ (10) 



Q j dp Wqpm{p) - Q j dp Wg; 



of integrals like in Eq. (10), described in Appendix A 

3. Results 

3.1. Statics 

It will be of further importance to understand the small-wave-number limits of the 
vertex Wqp. This depends on the dimensionality; for d = 3 one obtains from Eq. ([6]) 

I i_ n+p 

Considering p — )• at fixed g, {2iTYWqp = 2q'^v{q) + 0{p'^), and for g — at fixed p, 
{2'KYWqp = ^pv{p) +0{q'^). Thus, Wqp is regular in the g — >■ and p — >■ limits, and 
falls off quickly enough at large p to ensure the finiteness of the integrals we consider. 



Wqp^j—^j-^l kdk{q^ -p'^ + efv{k). (11) 



The infiuence of the small- wavenumber regime in Eq. ( 10 ) is easiest understood 



by introducing an infrared cutoff parameter 5 in the integral and discussing the 



asymptotic dependence on it. Close to the transition, m{q) is small; Eq. ( 10 1 becomes 
a linear integral equation with positive definite kernel Wqp. Nontrivial solutions only 
exist if Q is an eigenvalue |42J. The equivalent of the Perron- Frobenius theorem for 
nonnegative matrices holds, called Jentzsch's theorem in the context of Fredholm 
integral equations |43j : there exists a non-degenerate minimal eigenvalue Qc that 
is positive and real, and the corresponding eigenfunction can be chosen real and 
positive. Obviously, Qc defines the critical density: for densities below, only the trivial 
solution ms{q) — exists, so that the particle is delocalized. We denote the critical 
eigenfunction by h(q), 



hsiq) ^ Qc J dpWqphs{p). (12) 
We will also need the eigenfunction h{q) of the transposed equation, 

/■oo 

~hs{p) ^ Qc I dqhs{q)Wqp. (13) 



5 
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Figure 1. Critical eigenfunctions hs{q) (upper panel) and hj{q) (lower panel) 
for the Lorentz-gas MCT vertex for various S as indicated, normalized according 
to Eqs. | |14| l. The functions have been analytically extended to g < 5. The insets 
show hg {q) and hg (q) normalized to their respective maxima. 



In order to normalize the eigenfunctions we can impose two conditions, which we 
choose for later convenience, 

dq hs{q)hs{q) = 1 , hsmsiSf^^T^ ^ = 1 . (14) 



5 



Critical eigenfunctions hs{q) and hs{q) are shown in Fig. [l]for various cutoff values 
S, evaluated for the Lorentz-Gas vertex. Apart from the normalization, the shape of 
these functions depends rather insensitively on S, as exhibited by the insets of Fig. [l] 

The next-to-leading order expansion of Eg . (|10[) depends rather more sensitively 
on S. We shall distinguish "large" cutoff, S ^ y^m{S), from "small" ones, 6 <C y^'m{S). 



In the former case, the term m{p) in the denominator of the second part of Eq. ( 10 1 
is not relevant close to the transition, and the corresponding integral yields a term of 
C(m^). We thus expand for g> Qc 

m{q)^ahs{q)+a^Ks{q)+0{a'^), (5 large, (15) 

where cr is a small parameter quantifying the distance to the transition point in a way 



to be specified. Introducing e = {g — Qc)l Qc and inserting into Eq. (10 1 we get after 



multiplication with (q) and integration over q the solubility condition 

dqhs{q)ihsiq)/qy] ^ e 6 / {hs{S)hji5)) , 5 1arge,(16) 



s 



where in the approximation we have exploited that the integral is dominated by 
contributions from q 6. In other words, cr ~ e, which is the asymptotic relation also 
known from the schematic model |27] . 
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Figure 2. Order parameter m(0) as a function of density g, for various S as 
indicated. The inset shows the variation of the critical point, gc, with S. 



For small cutoff, 6 ^ y/ms{6), the second part in Eq. ( 10 1 needs to be treated with 
care. Specifically, in the small-p region we can replace Wqp « Wqs and ms{p) ~ ms{6) 
in order to perform this part of the integral explicitly. Omitting the remaining terms 
of O (ra^ ) , we arrive at 

/>oo 

m(p) = g dpWqpm{p) - gWgsm^^^S)'^ + 0{m^) . (17) 

We now need to expand differently. Again for g > Qc, 

m{q) = ahsiq) + a^''^ks{q) + ©(fJ^) , 5 small. (18) 

Inserting this into Eq. (17 1 and collecting terms of C'(ct'^/^), we obtain 



Ks{q) ^ gc j dpWqpKsip) (19) 
+ a-^'^eg, [ dpW^phsip) - gcWqshf\S)iTr/2) , (20) 



if we set e ~ a^^"^. The solubility condition is again obtained by contracting with 
hs{q) and reads 

^ = \e\/hs{S)hs{6f^\7T/2) , ,5 small, (21) 

where it is understood that e ^ corresponds to cr ^ and denotes the localized 
respective delocalized side of the transition. 

Thus, the asymptotic behavior of ms{q) changes depending on 6. Stated 
differently, for any fixed infrared cutoff there exist regimes close to the transition 
where the "large-(5" asymptote is reached. For distances e larger than a cross-over 
distance, the "small-(5" asymptote holds. We estimate this cross-over by equating the 
leading-order terms, from which one gets 



dq hs{q)hs{qY / q^ 



hs{5)hs{5){^/2Y5. (22) 
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To demonstrate this behavior, let us investigate the order parameter rn(0) as 
a function of density. Numerical results are shown in Fig. [2j Except in a region 
of densities around 0.6, the curves depend relatively weakly on S, indicating that the 
theory in d = 3 is for this aspect not sensitive to an infrared cutoff sufficiently far in the 



localized region. A numerical determination of the critical eigenvalue in Eq. ( 12 ) for 
the LG vertex gives Qc ~ 0.5644/i?'^ for 5 — )■ 0. As the inset of Fig. [2] demonstrates, for 
the d = 3 LG vertex, the variation of Qc with S is small. The localization point found in 
simulations of the three-dimensional Lorentz gas is Qc ~ 0.84/i?^ It is well known 

that MCT over-estimates the tendency to arrest, so that we will not be concerned 
by this discrepancy. Note that a generalized hydrodynamic approach, replacing all 

limit, allows to directly 
9/(47ri?3) 



occurrences of the memory kernel in Eq. ( 10 ) by its q 
evaluate the integral using Eq, 



_ , resuhing in = 9/(47ri?^) w 0.716R~^. 
The different asymptotic regimes are exhibited in Fig. [sj where rn(0) is shown as 
a function of e in double-logarithmic fashion. The asymptotes, Eq. ( |15[ ) together with 
Eq. (16), and Ec^. (18) with Eq. (21), are shown as dotted lines. They are found to 



agree well with the numerical results. Stars in Fig. [3] mark the cross-over point £*, 
Eq. (22), where the two asymptotic regimes meet. 

Figure |4] shows the noncrgodicity parameters f{q) obtained for S — 10~^ and 
6 — 0.01 for various e. Inserting the leading-order result m{q) = \a\hs{q) into Eq. ([9|, 
we get for a > 

/(9)-^^7^^^> (23) 



1 



qy{ahs{q))' 



so that for small q, f{q) becomes a Lorentzian from which we read off a length scale 
diverging as ^ a^^^^. 



3.2. Dynamical Scaling Functions 

We are now interested in the dynamic scaling of correlation functions in the limit of 
large times, i.e., time scales much longer than that of binary collisions. We shall thus 




Figure 4. Nonergodicity parameter f{q) corresponding to <5 = 10~* (solid lines) 
and 5 = 0.01 (dashed li nes) , for e = 10~"/^ as labeled. The asymptotic results 
given m(q)/q-^ and Eq. ||23| are shown in dotted and dash-dotted. 



analyze Eq. ^ in the limit 

zm{q, z) = 



\z\ V. Analogously to Eq. 

Q I dpWqpZ^l{p,z) + Q 



dpW, 



we get 

{zfi{p, z))^ 



IP 7 



(24) 



where ziJ,{p, z) = z'^ /v'^^+p'^iiyz + zm{p, z). 



s P- ~ zfi{p, z) 

Note in passing that the emerging integrals 
are well-defined if one assumes t) and 'm{q, t) to be completely monotone functions. 
In particular, there then holds for the real part dim{q, z) > for diz < j^, ensuring 
the analog of Eq. (17) with m replaced by fi on the r.h.s. One readily obtains that in 
leading order at the critical point, m{q, z) solves the eigenvalue equation ( 12 ) as 2; — 0. 
Thus, m{q,z) — hs{q)G{z) + 6m{q,z). The scaling ansatz observing m{q) ^ Odcrl) 
sets G{z) = \a\t^g{z) for rescaled frequencies z = zt^. Here is a time scale that 
diverges faster than l/\a\ as the transition is approached. The leading-order term 
is reminiscent of the standard MCT factorization theorem |21j : wave- number- and 
frequency-dependence factorize. For the correlator this implies (j>{q, t) ^ {h{q) / q^)G{t) 
as cr — > 0. But since t) < 1, at the localization transition the factorization theorem 
for 0((7, t) is violated even as t -> 00 for small q. This is demonstrated in Fig. |4j where 
dotted lines show f{q) ^ (jh{q)/q^ ~ m{q)/q^. 

To obtain the scaling equation, we need to consider the next-to-leading terms and 
hence distinguish again between small and large 5. In the latter regime, we set, in 
generalization of Eq. ( [T5| , 

Inserting into Eq. (24) and balancing the terms of 0(17^) 
the left-eigenfunction hs{qyi\ie scaling equation i{T /t„)z-\-£\a\zg{z)+e(j{zg{z))'^ = 0. 
Here we have used Eqs. (14) and (16) and defined t — v J dp hs{p)p^ . We thus fix 



K{q,z) + Oi\a\^), (5 large. (25) 
yields, after contraction with 



and get the two scaling solutions for cr ^ 0, 
= ±1 + g±iz) ± z{g±iz)f , 



S large. 



(26) 



(27) 
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Figure 5. Memory kernel m{q, t) at g = 5, for the case of large 5 (top panel) 
and small <5 (bottom panel), for distances to the respective transition points 
e = ±10~"/^ (ri = 1, . . . 8) and e = (curves labeled c). The asymptotic results 
obtained from the scaling equations, Eqs. | |27[ l and \31\ , axe shown as dash-dotted 
lines. Stars mark the time scales ta used in the scaling equation. Insets: m{q, t) 
on a linear-log scale. 



This is the scaUng equation derived earher by Gotze for the Fl model [H]. Its solution 
can be given as a closed expression for the inverse Laplace transform, g{t), where 
t = t/tcj. Here we are only interested in the asymptotes. Letting z — >■ 0, we find 
g+{z) ~ 1/z. Thus, g+{t) — > 1 for e > 0, i.e., the signature of localized dynamics. On 
the other hand g-{z) remains regular, indicating that g-{t) decays faster than 1/t at 
long times (exponential in fact). Letting 5 — >■ oo we obtain g{t) ~ (tt^)"^/^, so that 
for the time window l/v t <^ to-, we have the cr-independent asymptote 

-(.,^)-H'#(^)"''-^47^ ^ large. (28) 

Note that an explicit dependence on the cutoff parameter remains in the prefactor. 
The top panel of Fig. [sjexemplifies m{q = 5, t) in the large-(5 regime for various distances 
£ to the transition. Dash-dotted lines indicate the above asymptotes. 

Along similar lines, we deduce the small-5 scaling function. Here we have to 
expand as in Eq. (18 1, 

m{q, z) = \a\tcrhs{q)g{z) -f \af^'^t„k{q, z) , b small. (29) 

Balancing the lup^^ terms in Eq. (24) we obtain after multiplying with /15(g) and 
integration over q 

i{T/t,)z + e\a\zg{z) + h,{5)hs(5f'\^ l2)\a\^'Hzg{z)f'^ = . (30) 



Now, using Eqs. \2\\ and ( [26^ , this yields for e ^ 

Q = z{g±{z)f + {g±{z)f±2ig±{z)-l, J small (31) 

As i — > 00, this equation is solved by g±{z) = i{iz)^^/^ , and thus g± (i) ^ {^^/■^/r(l/3) 
for l/i' <^ t <^ ta- As £ — >■ 0, we note again g+{z) ~ —'i-fz (locahzed solution), but 
g-{z) ^ i — leading to g~-(i) ~ / y/A^. The behavior we find is thus 

m(g,i) - (/i5(5)7r/2)-2/3/r(i/3)(t/r)-2/3^ f < t^, 5 small, (32) 

m{q,t) -^\a\hs{q)/y^{t/tay^/'^ , t > i,, , e ^- 0" , (5 small (33) 
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Figure 6. m{5, t) at the critical point as a function oi t/t* according to Eq. | |34| l 
for various <5 as indicated. For the unlabeled curve, 5 = 10~^/R. 



These power laws are exemplified in the lower panel of Fig. [sj for various e in the small- 
6 regime. As for the large-5 regime (upper panel), the asymptotic results describe the 
numerical solution for m{q,t) well even for e = C(O.l). Treating the small-wave- 
number regime of the vertex seriously, a long-time tail for m{q, t) results; however, 
this is not the known universal long-time tail of the LG [3HI3S] in d = 3, t~^/^, which 
is not contained in the present MCT. 

The two regimes outlined above, large and small S are separated by the condition 
6 ~ ^y\zm{S, z)\. For fixed S, this implies a time scale t* where the small-(5 asymptotics 
for t <^t* crosses over to the Xaige-S asymptotics for t ^ t* . Equating Eqs. (28) and 



(32 1, and using the approximation Eq. (22), we find 

16 T 



TTT{l/3fhs{S) 



(34) 



Thus, plotting ■m{6,t) at various 5 as a function of t ~ t/t* reveals a master curve 
with asymptotes i^^/^ for t < 1, and for i > 1, for cr 0; see Fig. [i] 

If the cutoff is held fixed, the time scale io- marks the cross-over between the 
different power-law regimes. The small-J asymptotes are always encountered first, and 
if the corresponding t^ < t* , the cutoff- influenced large-5 asymptote is not seen. Only 
close enough to the transition does the cutoff-dependent small-(5 asymptote become 
evident for long times. This can be seen in Fig.[7j where a star marks t* , and dashed 
and dash-dotted lines the different asymptotic solutions. 



3.3. Mean- Squared Displacement 

The mean-squared displacement 5r'^{t) is obtained from the q — >■ limit of the 
correlator, 4>{q,t) = 1 — {q^ /6)6r^{t). Inserting into Eq. ([sj, one gets after Laplace- 
transforming the equivalent of Eq. ([2| , 

^' z + iiy + m{0, z) \a\hs{0)t^ Pg{z) ' ^ ' 

The latter asymptotic behavior is seen from inserting the leading-order result for 
m(0, z) and noticing that io- — > oo close to the transition. The generic features of 
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Figure 7. m{S,t) as a function of t, for 5 = 10~^, at various distances to the 
transition, £ = ±10""/^ as labeled and e = (label c). Dashed lines show small-(5 
asymptotes, dash-dotted lines the large-<5 asymptotes at the critical point and for 
exemplary e. A star marks t* , Eq. ||34|l. 



Sr'^{t) at long times are then immediately clear from those of m(0,t). For e > 0, from 
g{z) ^ — l/z there results the localization length 



= \/ Sr^{oo)/6 = 



\a\^/^hy\0) 



e>0. 



(36) 



On the liquid side, g-{z) remains regular as z 

t 



Sr''{t oo) 



■ 0. We then get ordinary diffusion, 
e < 0. (37) 



\a\hs{0) 

In the intermediate time window, X/v t <^ t„, the power laws identified for the 
memory kernel, m{q^t) ^ with < a; < 1, result in power laws (5r^(t) ^ P. We get 

(38) 



(5r2(t) 



an intermediate window of anomalous power-law diffusion, with exponents a; = 1/2 in 
the large-5 case, and a; = 2/3 in the case of small 5. 

Figure [s] exemplifies the behavior for 5 = 10~^, corresponding to the memory 
kernels shown in Fig. [7j As stated above, the results are for overdamped Brownian 
dynamics; this version of Eq. (35) reads (5r^(z) = (6/z^)/[w-|-m(0, z)], and the relevant 
asymptotic results are obtained by replacing v"^^ = 1 in Eqs. ( 36 1 to ( 38 ) . Again, a time 
scale t* can be identified where the contributions from the two asymptotic solutions 
balance. As in Eq. (34), t*/r ^ l/tJ"^, but the prefactors are different: 

6 



t* 



2 / 2V2r(5/3) 



r 

(53 



(39) 



Note that t* sa lOt*, as is verified by comparing the marked positions in Fig. [t] and [s] 
Recall that depending on 5 the connection between e and a changes. However, 
\a\tcr — T/\e\, and ^.^(O) changes relatively weakly with 5, so that Eq. ( [ST] ) yields the 
same diffusion law, 6r'^{t) = 6Dt irrespective of S, with a diffusion coefficient D ~ \e\. 
On the localized side of the transition, Eq. (36) yields Tc ~ for the small-5 case, 
but Tc ~ e~^/^ for the large-5 case. 
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4. Conclusions 

We have discussed the different asymptotic regimes that emerge from MCT for the 
tagged-particle dynamics close to a localization transition, depending on an arbitrarily 
chosen infrared cutoff parameter S. Such a parameter is not justified a priori, but 
can be rationalized by recognizing that the small-g structure of the MCT vertex for 
the LG case studied here is erroneous [28|. Reassuringly, in d = 3, the predicted 
position of the localization transition depends only weakly on S, although its overall 
value suffers from similar disagreement with computer-simulation result as the MCT- 
predicted glass transition point. Note however that in c? = 2 the situation changes. 
Following an argument of Ref. 46J , one can show that the present MCT does not allow 
solutions recovering long-time diffusion if (5 — >■ 0; this results in a transition point that 
shifts to zero with decreasing S. An analysis like the one presented here can still be 
carried through, if the changed structure of the vertex Vqp is taken into account |47j . 

Depending on the distance to the transition e, and the cutoff S, one observes 
power laws with different exponents for the memory kernel m{q,t); one distinguishes 
the small-(5 dynamics from the large-(5 case, separated by a time scale t* . This also 
holds for other dynamic quantities, like the mean-squared displacement Sr'^{t). The 
latter grows as t^^^ for the large-;? case, which is the result also known from schematic 
MCT (27j and is the relevant result for previous numerical MCT calculations that do 
not treat the small-g anomalies of the vertex seriously. It is also, together with the 
D ~ \e\ obtained above, the result expected from mean-field theory f48^. The results 
presented here also serve the purpose to highlight which parameter regimes can safely 
be treated by the present MCT in its widely used incarnation on a discrete wave- vector 
grid (where typically, 6 « 0.2 /R), at least in c? = 3. 

Our introduction of a cutoff has analogies to a modern formulation of the 
renormalization group (RG) approach where an infra-red cutoff controls the RG 
flow of the effective average action. It would be promising to establish an RG flow for 
the MCT equations, possibly following along the lines of a cluster-MCT ^50j, although 
it is not obvious how to do this. 



REFERENCES 



14 



Acknowledgment s 

This work was funded in part by the Deutsche Forschungsgemeinschaft, DFG, project 
P8 of the Research Unit FOR1394 "NonUnear Response to Probe Vitrification". 
Th. V. thanks for funding through the Hchxiholtz-Forschungsgemeinschaft (Impuls- 
und Vernetzungsfonds, VH-NG-406) and ZukunftskoUeg, Universitat Konstanz. 



Appendix A. Numerical Procedure 



We sketch the numerical scheme used to solve Eq. ( 10 ) and similar integral equations. 
To capture the behavior at small wave numbers, a logarithmic grid was introduced, 
containing Ns points starting at S up to some Aq; from Aq to an upper cutoff Q, 
Nq equidistantly spaced points complete the grid. We used Ns = 50, Nq = 300, 
Aq = 0.08/i?, and Q = 48/i? for static calcuations. Time-dependent quantities were 
obtained with Aq — 0.4 and Q = 2A/R. These discretizations were found to be 
sufficiently accurate and have been checked with larger grids in some cases [17]. 

The critical eigenfunctions hs[q) and hs[q) have been obtained by iteratively 



solving Eq. ( 10 ) starting with some constant initial guess; at each iteration step the 
approximate solution is inserted in the r.h.s., and the result is taken as the next 
approximant. As the iteration convergence is dominated by the eigenvalue Qc, this 
procedure also determines the critical density, calculated as the ratio between the 
results of two iteration steps at fixed q. 

For the time-dependent Eq. ([s]), a standard MCT algorithm is used where an 
equidistant grid is chosen for the time i, which is repeatedly coarsened by a factor 2 
after a fixed number of steps in order to cover the required number of decades. 
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